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Abstract 

We consider an Einstein-Hilbert-Dilaton action for gravity coupled to var- 
ious types of Abelian and non-Abelian gauge fields in a spatially finite sys- 
tem. These include Yang-Mills fields and Abelian gauge fields with three 
and four-form field strengths. We obtain various quasilocal quantities as- 
sociated with these fields, including their energy and angular momentum, 
and develop methods for calculating conserved charges when a solution 
possesses sufficient symmetry. For stationary black holes, we find an ex- 
pression for the entropy from the micro-canonical form of the action. We 
also find a form of the first law of black hole thermodynamics for black 
holes with the gauge fields of the type considered here. 
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1 Introduction 



The relationship between the Euclidean-action formulation of quantum gravity 
and the thermodynamics of the gravitational field has been a subject of increas- 
ing interest in recent years. Fundamental connections between the partition 
function of the grand canonical ensemble and the Euclidean-action path inte- 
gral were first pointed out by Gibbons and Hawking Q], who argued that the 
Euclidean gravitational action is equal to the grand canonical free energy times 
the reciprocal of the temperature associated with a black hole (or cosmological) 
event horizon 0]. 

More recently Brown and York have extended this work by considering the 
formulation of the partition function for gravitating systems of finite spatial 
extent j^, ||. Virtually all systems with which we have any experience have a 
finite spatial boundary; indeed, one of the central concepts in thermodynamics 
is that of a system and a reservoir that are separated by a partition. A phys- 
ical realization oLihese concepts is needed in order to apply thermodynamics 
in a sensible way.u When we use thermodynamics to describe self-gravitating 
systems (including objects such as black holes), we must divide our space-time 
into a region which contains the system (the black hole) and the remainder of 
the space-time which can be treated as the reservoir. The traditional practice is 
to study black hole thermodynamics at space-like infinity, assuming reasonable 
asymptotic conditions without any boundary. However, this approach has a 
number of deficiencies. First, it requires that a spacetime which is a solution 
to the field equations (and often the accompanying matter matter fields) pos- 
sess appropriate asymptotic behaviour, typically asymptotic flatness. However 
asymptotic flatness is never satisfied in reality, and is not always an appropriate 
theoretical idealization: many black hole spacetimes exist that are solutions to 
the (dilatonic) gravitational field equations that do not possess asymptotic flat- 
ness. Some have been found recently that are not even asymptotic to de Sitter 
or anti-de Sitter spacetime 0. Furthermore, in the study of pair creation of 
black holes one is forced to consider non-asymptotically flat spacetimes with an 
acceleration horizon [[}] , necessitating a more careful consideration of boundary 
terms in the formulation of the Hamiltonian 0. Second, construction of a par- 
tition function (which is central in the study of statistical mechanics) requires 
the stability of the system, which is only realized when a finite size is imposed 
[0, a point also noted by Hay ward and Wong B. For example, the heat ca- 
pacity for a Schwarzschild black hole is negative B] if one fixes the temperature 
at infinity, and the formal expression for the partition function is not logically 



consistent 10 . However if the temperature is fixed at a finite spatial boundary, 



1 Another central concept of thermodynamics is, of course, thermodynamic equilibrium. 
The realization of this is the stationarity of a system, where we say a system is stationary 
if there exists a time-like (Killing) vector field for which the Lie derivative of all the fields 
considered vanish, and the boundary of the system is chosen to contain the orbits of these 
vectors. We shall assume that our systems are stationary. 
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there is no inconsistency in the black hole partition function and the heat capac- 
ity is positive O] . This approach may also be extended to black holes in anti 
de Sitter spacetimes 12 , where analogous results are obtained. Finally, it seems 
to us that, on physical grounds, one should be able to define thermodynamics 
appropriate to observers who are at a finite distance from the black hole. 

It is important, therefore, to construct thermodynamic quantities appropri- 
ate to observers at the (finite) boundary of a system. The quasilocal formalism, 
which has been developed extensively by Brown and York [fl3| provides us with 
a means. This formalism is based on a Hamilton-Jacobi principle wherein the 
boundary terms of the action functional for a compact region give rise to the 
quasilocal quantities such as energy and angular momentum. Brown and York 
[|| have also developed a paradigm for understanding the relationship between 
the classical mechanics, the statistical mechanics, and the thermodynamics of 
gravitating systems, based on the boundary conditions of the system. Here, the 
fields that are held fixed on the boundary in deriving the classical mechanics 
equations of motion from an action principle determine the statistical ensemble 
of the corresponding statistical mechanics, and thus the type of thermodynamic 
partition that must be imposed. There is a parallelism between the Legendre 
transformation — a canonical transformation — of the boundary terms of the ac- 
tion, and the Laplace transformation — which changes the type of ensemble — of 
the path integral for the statistical mechanics. This parallelism may be used 
to identify a "micro-canonical" action that has boundary terms appropriate to 
a micro-canonical statistical ensemble ||, and from it an expression for the 
entropy of a thermodynamic system may be obtained. 

In this paper we extend this formalism to include the most general action of 
gauge fields coupled to dilaton gravity that has at most two derivatives in any 
term. This class of actions includes the low-energy limit to string theories [|14| 
and is also of interest as an empirical foil for testing general relativity JIB], |l(i| . 
We consider non-Abelian gauge theories coupled to gravity, the dilaton, and an 
axion field, as well as Abelian two-form and three-form gauge potentials with 
three-form and four-form field strengths respectively, also coupled to a dilation. 
(We note that recent research in non-Abelian gauge fields coupled to gravity has 
led to new black hole solutions that refute much of the folklore about black holes, 
such as the "no hair" conjecture [0). We will find that a form of the "first law" 
of black hole thermodynamics can be reconstructed for the theory considered 
which is quite reminiscent of the conventional one for the Einstein-Maxwell 
theory. Although we restrict ourselves to four dimensional space-times, higher 
dimensional generalizations of our work are straightforward. For each sector of 
the action, we obtain expressions for the quasilocal thermodynamic quantities 
that will appear in the first law of thermodynamics. In addition, we discuss the 
role of conserved quantities associated with the gauge and gravitational fields. In 
general, the presence of a conserved quantity depends on the type of solution. 
Such quantities always require the presence of some sort of symmetry in the 
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solution. (Solutions that satisfy certain asymptotic conditions often also possess 
an asymptotic symmetry that can be used to define conserved quantities. Here, 
the solution must at least possess exact symmetries on some finite boundary in 
order to admit conserved quantities for finite-sized systems.) We also give a brief 
summary of the construction of a statistical mechanics for the class of theories 
considered, and we obtain an expression for the entropy and a form of the first 
law of thermodynamics for systems possessing an event horizon. We note that a 
version of the first law of black hole thermodynamics for a finite system has been 
obtained for vacuum general relativity as well as vacuum Einstein-Maxwell fl8f ; 
we recover these results in the case of an Abelian gauge theory with vanishing 
dilaton and axion fields. 

Closely related to the work of Brown and York, as well as to our present work, 
is that of Wald's Noether charge formalism (l9). Wald obtains an expression 
for the entropy of a space-time solution of a very general class of theories that 
are required to be derivable via an action principle from a Lagrangian density 
that is covariant under diffcomorphisms. Here, the entropy is identified as the 
Noether charge associated with this covariance. Recently, Iyer and Wald [^0| 
have extended their treatment to systems with finite boundaries, and they have 
shown how to formally recover some of the results of Brown and York |l3| . The 
present work can be viewed, therefore, as complementary to this technique. We 
restrict our considerations to actions with at most two derivatives in every term, 
and we explicity evaluate quantities which are only implicitly defined in ref. 
for specific types of boundary conditions. We justify the choices of boundary 
conditions by appealing to the role they play in the connection of the classical 
mechanics, the statistical mechanics, and the thermodynamics discussed above. 

Therefore, our agenda is the following: In section two, we will study the 
Einstein-Hilbert-Dilaton sector, which we will consider as our gravitational the- 
ory in the absence of additional fields. In section three, we consider a Yang-Mills 
field that couples to both the metric and the dilaton, extending this to include 
axion couplings in section four. In section five, we turn to an Abelian gauge the- 
ory involving a four-form field strength (also coupled to the dilaton), and briefly 
consider its relation to a cosmological constant. In section six we complete our 
survey of matter fields with an Abelian three-form field strength gauge theory 
coupled to the dilation. Section seven is a review of statistical mechanics based 
on path integral techniques. Here we adopt a "general form" for the gauge fields 
rather than treating each case separately. In section eight we apply these results 
to obtain an expression for the entropy of a system containing a black hole and a 
form of the first law of black hole thermodynamics. Concluding remarks follow. 

In this paper, we adopt the conventions of Wald |^l[ . Units are chosen so 
that the speed of light, Newton's constant and the rationalized Planck constant 
are all taken to be unity. A summary of the notation of the manifolds, fields, 
and related quantities considered in this paper is given in the appendix. 
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2 Einstein-Hilbert-Dilaton Sector 

The theory of gravity that we consider is one which is based on the usual action 
of general relativity, but with the addition of a scalar function that couples to 
the curvature called the dilaton. In general, it is possible to redefine the metric 
via a conformal transformation so that the dilaton does not couple to the new 
curvature. However, as we shall see, such a transformation affects physical 
quantities such as the entropy associated with black holes, so we consider the 
more general case. 

As indicated, our theory of gravitation will be based on an action principle, 
and is closely related to the usual Einstcin-Hilbert action. The gravitational 
field equations may be deduced using Hamilton's principle when considering 
variations in the geometry. It is now known that serious restrictions must be 
placed on the boundary of the space-time region in order that such a variational 
principle be well defined. In particular, both the variation of the induced metric 
and its derivatives must be held fixed on the boundary. Alternately, the action 
may be supplemented with additional boundary terms such that we need only 
fix the induced metric on the boundary. We take the latter approach; from these 
boundary terms many useful quantities can be defined, as we will show. 

In this section, we look at the gravitational sector of the theory, and we 
defer consideration of various types of matter to later sections. Insofar as we 
restrict ourselves to dilatonic gravity, we will consider a very general action. The 
coupling to the curvature, the kinetic energy of the dilaton, and the potential 
energy of the dilaton will all be arbitrary functions of the dilaton (alone) . This 
practice will be continued in the next section when we consider gauge fields with 
arbitrary couplings to the dilaton. 

We will restrict our considerations to a spatially finite region of a four- 
dimensional manifold, M. , which has a topology of £ x H where £ is a space-like 
hypersurface and IR is a real interval. On this manifold, we define a metric, g^, 
and its compatible derivative operator V M . Objects with Greek indices represent 
tensor quantities on the four dimensional manifold. We will consider two parts 
of the boundary of this manifold: the outer boundary 7^ ut0 r = d£ x I, and the 
initial space- like hypersurface initial- The addition of an inner boundary and 
a final space-like hypersurface are trivial extensions of the present analysis that 
will be important in following sections. Here we shall refer to them simply as 
T and £ respectively. 

Associated with T is an outward directed space-like normal vector n^. We 
can construct the first and second fundamental forms 7^ = g^ — n^riy and 
O^j = —\£ n lnv These are to be viewed as tensors on T, and the indices i, j, 
etc., will denote tensor quantities on T. However, we will use j?, defined in an 
analogous manner to 7^, to denote the projection operator from M. onto T. 
The derivative operator compatible with 7^ is Aj. 

We can foliate M. into space-like leaves £ t via the parameter t which is a 
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coordinate along R; associate with it the vector field t M = d^t. View the initial 
hypersurface as one of these leaves, say, the t — one. Little confusion arises 
from dropping the t index and considering an arbitrary leaf, S, in the foliation. 
There is a future directed time-like normal vector to S. We may define the 
fundamental forms, — g^ + u^u v and K^ u — —\£, u h ilv — these are viewed 
as tensors on S, and such tensors will be shown with indices i, j, etc. — and 
projection operator, from M onto S. We may also define the lapse and 
shift of the foliation by N = -t^u^ and N l = ti^, so that f = Nu^ + JV>. 
The derivative operator on S compatible with is given by V», and this can 
be used to define the Ricci scalar on S: R[h]. It is also important to consider the 
boundary, dS of S. We require that, on the intersection of S and T, u^n^ = 0. 
The first and second fundamental forms on OS are <Tjj = — Hfiij and fcjj = 
— ^£ n (jjj. The Latin indices a, b, etc., are used to denote tensors on dS. The 
projection operator from T onto dS, <j l a , is obtained from = jij + UiUj. 

A zero vorticity observer is one for whom the vorticity, vj = *(v A dv), 
is zero, where v is the velocity of the observer. It can be seen that observers 
who are co-moving with a given foliation on the boundary T are zero vorticity 
observers. For such an observer, the velocity is just the normal vector u M , and 
the acceleration of the normal vector is given by = u u V v u^ = N~ 1 WV V N . 

A summary of the notation described above is given in Table 1 of the ap- 
pendix. 

2.1 Variation of the Action 

In accordance with the above considerations, we choose the action for the grav- 
itational sector to be: 

"5*EHD = 

Jm 

-2 J^d 3 x^f EK (#) tr(0) 

-2 f d 3 xVhf EH (V)tr(K) (2.1) 

where the dilaton field is given by •f', and / EH (!f - ), / KE (^), and / PE (^) are model- 
dependent functions of the dilaton. They are arbitrary insofar as they contain 
no fields other than the dilaton and no derivatives of the dilaton. Notice the 
presence of the two boundary terms: these are just what is needed to make a 
well defined variational principle on the initial space-like hypersurface and outer 
boundary as we shall see below. 

Now we vary the geometry and dilaton field configurations. The geometry 
is varied subject to the (gauge) restriction that the leaves of foliation remain 
orthogonal to the boundary. (That is, we hold the boundary fixed in that 
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variations of the normal dual- vectors to the boundaries are proportional to the 
normal dual- vectors.) For convenience, we vary the inverse metric, except on the 
boundaries where we can write the variation in terms of the covariant induced 
metrics. The induced variation in the action is given by: 

J M 

d 3 x (7r«*y« + /ZdiiMO 

T 

d 3 x (p^Shij + P A a6W) . (2.2) 

E 



Here, 



-2(V^)(V^)) + |^/ M (!P) (2.3) 
(with G/tulg] = R^ulg] - ^g^uRlg]) and 

Sdil = + - 2Vl/ n (W) + ^ (2.4) 

can be considered as equations in the following way: under variations of the 
geometry and dilaton field that leave the geometry and dilaton configurations 
on the boundary fixed, (Sgcom)^!/ and S^n will be zero at an extremum of the 
action. (In general, there will be contributions to these equations of motion 
from any additional matter present.) 

Alternately, we could consider variations of the geometry and dilaton con- 
figurations in which the equations of motion are held fixed, but the boundary 
geometry and dilaton configurations are varied. Under these variations, we see 



from equation (2.2) that we can define the momenta conjugate to the boundary 
geometry of the initial and outer hypersurfaces as 7T lJ = (^Sehd/^TijOcl and 
= —(5S EHD /5hij) clj respectively. (The subscripted 'cl' emphasizes that the 
variations are amongst field configurations satisfying the equations of motion.) 
Explicitly, these are: 

n ij = V^l (f 3 '(n^Um) + IU*){® lJ - l 13 tr(0))) (2.5) 



and 



-Vh (^(^/ M (ff)) + M^K* - h"ix(K))) (2.6) 

Similarly, there are momenta conjugate to the dilaton field on the initial and 
outer hypersurfaces defined in an analogous way. These quantities are: 

n ™ = (^PL = ^ (2/-(*K»M* - tr(0)) (2.7) 
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and 

ftn = -(^) CL - -^(2/ KE («^2^tr(if)). (2.8) 
2.2 Quasilocal Quantities 

The T boundary momentum, tv ij , contains useful information about the energy 
and momentum densities of the gravitational field held within this boundary. To 
identify these quantities, it is useful to decompose the variation of the boundary 
metric, jij, into the various projections normal and onto the foliation. Thus, 
we write: 

&Hi = dftfSaoi, - —8Nu lU] - —u {l a j)a 5N a . (2.9) 

The corresponding decomposition of the boundary momentum, Tt lJ , conjugate to 
7ij leads us to define the surface energy density, the surface momentum density, 
and the surface stress density as: 

£ = ^UiUjT ij , (Jehd) =-V^u^t 11 , and S ab = y/ffafa]^ 

respectively, where t ij — 2it lJ / 7 is the surface stress-energy- momentum on 
T. 



We next use equation (2.5), and the relationship: 

©ij = kij + UiUj^a^) + u^cr^mPK^ ( 2 -10) 

(recall that a M = u u V „u M is the acceleration of the normal u 11 ) to obtain expres- 
sions for £, (Jehd)", and S ab . An immediate consequence of equation ( |2.10 ) is 
that tr(0) = tr(fc) — n^a^. Then, we have: 

£ = -2V^(n^^/ EH (^) - / EH (^)tr(fc)) - £ (2.11) 

(Jbhd)" = 2V^hn(&)n- % a?KV = -h^-^p* - (J ESD , ) a (2.12) 

and 

S ab = 2^d (a ab n^f w (&) 

+f EH (<P)(k ab -a ab (tr(k) -n^))) - (S„) ab (2.13) 

where £ , (iTehd.o) , and (S ) are additional contributions that arise in sup- 
plementing the action S E n-D with an additional (reference) functional of the 
boundary fields, as we will discuss in the following sub-section. 
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We can construct a quantity, y, from the momentum TTdii conjugate to 
the dilaton configuration on T: y — N^IIm — y o where y o is an arbitrary 
background contribution that will be discussed below. Thus we have: 

y = ^ (2/ KB (<^Kc^ - (tr(fc) - n"o p )) - y, (2.14) 

Note that y is also a sca lar de nsity on the two-surface dS. Using the definitions 
in equations (P3T| - ^13|) and ( |l4[) , we can write the T boundary terms in the 



variation of the action as: 

SS EHD \ T = I d 3 x (-SSN + (J EHU ) a SN a + N{\S ab 5a ab + ySV)). (2.15) 



It is worth noting that £, (J7ehd)ci, &ab, and W are all extensive variables 
because they can be constructed out of the phase space variables, (p' J , hjj) 
and (Pdii,^) on S. However the lapse and shift cannot be constructed out 
of this information, and such quantities are called intensive variables. We see 



that, in equation (2.15), the first two terms of the integrand involve variations 
of intensive variables (with extensive variables as coefficients), while the last 
three terms involve variations of extensive variables (with coefficients that are 
intensive variables due to the lapse function) . We will consider the implications 
of this below. 

Define the quasilocal energy to be the integral over the two-surface dS of 
the quasilocal energy density: 



E = / d z x£. (2.16) 

This quantity has useful properties such as additivity although it is not 
necessarily positive definite. The quasilocal energy is observer-dependent, i.e., 
even if we are given a natural choice of boundary T, the above definition of 
the quasilocal energy will still depend on the foliation; in particular, how the 
leaves £ intersect with T. However, we would expect that we should be able 
to do better than this when the space-time is stationary. The construction of 
conserved charges resulting from space-times with symmetry will be addressed 
below. 



2.3 Reference Space-time 

The complete gravitational action will involve the action S^hd given in equa- 



tion (2.1), plus some additional functional of the boundary metric and dilaton 
configuration. This additional piece, S , does not contribute to the equations 
of motion because it contributes only to the boundary. It does, however, con- 
tribute to the momenta conjugate to the gravitational and dilaton fields and 
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thus to the quasilocal quantities discussed above. It will be sufficient in the 
present analysis to consider S to be a functional on the boundary T alone. 
Furthermore, we will assume that it is a functional of the metric 7^, and not 
T boundary dilaton configuration, but we will not consider the general case in 
which it is also a functional of boundary configurations of matter fields (though 
the inclusion of these is straight-forward). 

The specification of the functional S is akin to the specification of some 
reference space-time. Although there are restrictions on the types of reference 
space-time allowed (which have to do with the embedding of the two surfaces 
dS in the reference spacetime), we will assume that a natural reference space- 
time can be found. Often, this is can be achieved simply by setting constants of 
integration of a particular solution to some special value that then specifies the 
reference. Our restriction that the functional S not be a functional of matter 
fields (other than the dilaton) is just a specification of a vacuum reference space- 
time. 

The primary advantage of the above interpretation of S is that it ensures 
that the action is a linear functional of lapse and shift and, as we shall see, this 
means that the extensive variables defined above continue to be extensive. (This 
is the principal restriction that must be placed on S .) We shall thus write the 
additional piece in the suggestive form: 

S = J d 3 x (N£ - N a (J EHD , ) a ). (2.17) 

Here, £ and ( l 7 E HD,o)a can be constructed out of the phase-space variables alone 
(and are thus extensive), and can be identified by the functional derivative of 
S with respect to the lapse and the (negative) shift respectively. The surface- 
stress and dilaton density of the reference spacetime can be obtained from the 
following relation: 

d 3 xN(±(S„) ab 6a ab +y„5V) = - [ d 3 x (NS£ Q - N a 5(J EliD , Q ) a ) (2.18) 

which allows one to determine (S ) ab and y o when the form of £ and ( l 7 EH D.o)a 
are known in terms of a a b and \P. The net result of the add ition of S to S Erm 
is the inclusion of the extra terms in equations ( 2.11 - 2~13| ) as well as in the 



definition of y. Of course, one possible choice is always S = in which case 
these terms would be absent. 



2.4 Conserved Charges 

Given a set of observers, all having histories on T, we wish to identify quantities 
related to the geometry of M. that are the same for all the observers. We will 
require the presence of a Killing vector, £*, on T. Our first step is to consider 
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the equation of motion for the geometry. We have been deferring a detailed 
discussion of matter to a later section; however, here we will allow for the 
presence of matter so that the equation of motion reads 

{~ > gcom)fiu 2^/41/ (2.19) 

where Tn V is the stress-energy-momentum tensor of the matter and (Hgcom)^ is 



given by equation (2.3). Computing n M 7^(^ gco -m)nv with the aid of the Gauss- 
Codacci relationship 

n»f»R pv [g] = - A, (0« - 7 « tr(6>)) (2.20) 

we obtain 

v^AjT* = iT dil A^ - n^^Tftv (2.21) 

which shows that the surface stress-energy-momentum tensor is not divergence- 
less: it has source terms arising from the presence of the dilaton and the matter. 



However, we can contract equation (2.21) with the Killing vector. In addition 



to the usual requirements of a Killing vector, we require that the dilaton is 
constant on orbits of the Killing vector. That is, we assume that satisfies 
both the usual Killing equation, £^jij — 0, for a Killing vector on T, as well as 
£^ = 0. The left hand side of equation ( [2.21 ) becomes a total divergence due 



to the Killing equation and the symmetry of t vj . Integrating over T, we find: 

'Tnj; in j tia , .'7 



/•Tni^final P 

/ d 2 xV^^iUjT ij = / (Pxsf^CTi^. (2.22) 



I n th e event that the right hand side vanishes for arbitrary -Efi na i, equation 
( [2.22 ) expresses a conservation law for a geometric charge, 



/ dir x n . (2.23) 

Jas 



ids 

Contingent upon the type of matter present, there are many reasons why this 
may be the case. First, it may be that T is positioned such that there is little 
matter in its vicinity, so is negligible. Secondly, it may be the case that 
the particular projection t^Ti^n^ vanishes. Alternately, if T^ v is conserved, 
that is V \ i T^ 1 ' — 0, and if the Killing vector field on T can be promoted to a 
Killing field over A4, (that is, if £ M is a Killing vector field on A4, satisfying the 
additional requirements above, and T contains the orbits of these vectors) then 
we can re-write equation ( |2.22j ) in the form: 



k(T n i: ini ti a - k(T n £ final ) = f fmal d 3 xVhu"CT^ (2.24) 

•^initial 



and the integrand on the right hand side of this equation may vanish. 
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Suppose if* is a space-like azimuthal Killing vector. Then, we can define an 
angular momentum as J = k[<p]. If the surface dS is taken so that it contains 
the orbits of the Killing vector, then we can write: 

J= / d 2 x(J EH ») a ip a . (2.25) 
J as 

When is time-like, we can define a mass as M = — k[£]. If the spacctimc 
is also static, that is, is surface forming, then we can choose a two-surface 
dS for which the Killing vector is proportional to the time-like normal. In this 
case, the mass can be written in the form: 

M= / d 2 xN£. (2.26) 

JdS 

Note that, in general, the quasilocal energy will not agree with the conserved 
mass of the space-time unless u % is a Killing vector on T. However, in asymptot- 
ically flat, static space-times for which T is taken at spatial infinity, a foliation 
in which u 1 approaches the time-like Killing vector is usually adopted, and then 
the definitions of mass and quasilocal energy given here will agree. 



2.5 Canonical Form of the Action 

We now turn to the canonical decomposition of the action of equation ( |2 . 1[ ) . 
The Hamiltonian density of the Einstein-Hilbert-Dilaton sector is defined as: 

H EnD = P ^£ t h,j + P dn £ t & - L EHD (2.27) 

where L EHD is the Lagrangian density which, when integrated over M. yields 
the action of equation ( 2A_ ) . To evaluate the first term, we use the definition of 
the second fundamental form, and the relationship between and f\ We find 
that: 

//'l',//., = -2iVp^ify - 2iV,Vj/ J + 2Vj(iV^"). (2.28) 
Similarly, the second term is 

P da £t$ - NP dil $ + N\P dil V^) (2.29) 

where the symbol >P is shorthand for u^d^. The Lagrangian density can be 
canonically decomposed using the Gauss-Codacci relationship: 

R[g] =R[h}+ KVKi, - (tr(Jf )) 2 - 2V P («" tr(K) + a") (2.30) 
as well as 

(W) 2 = (V^) 2 -4> 2 . (2.31) 
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We also use the relationship = N 1 VjA r which is appropriate for zero- 
vorticity observers. The Hamiltonian can then be written as: 

d x -/? EHD 

s 



d 3 x (H EH »N + (n EUD )jN~ l ) 



a" 

E 

+ [ d 2 x(£N-(J EHD ) a N a ) (2.32) 

JdS 

where, 

H EHD = -2p^ + P dn !P 

-Vh (j EH {V )(R[h] + K^K- l3 - (tr(X)) : 

+2tr(X)^M^-2V 2 / EH (f) 
aw 

+/ KE (^)(V<^) 2 - f KE (&)& 2 + / PE (30) (2.33) 

is the Hamiltonian constraint and 

(^ehd)i = -2Vj^ + PdiiV^ (2.34) 

is the momentum constraint. When the constraint equations hold, the Hamil- 
tonian is purely a boundary term. 

The action, equation (2.1), can now be re-written in canonical form: 



$EHD — 



/ d 4 x( P ^£ t h, 3 + P dil £ t ^-H EHD ) 

J M 

= J dtf^J d 3 x ( P ^£ t h- l3 + P dU £ 4 «f - H EHD N - (WehdJ^/V) 

+ / d 2 x(-£N+(J EHD ) a N a )). (2.35) 

JdS / 

This form of the action will be useful later in the study of the thermodynamics 
of gravitational systems. 



3 Yang-Mills Sector 

Here we consider non-Abelian gauge fields with an arbitrary gauge group (S. 
The case of electro-magnetism is, of course, just a simplification of the general 
results of this section when the gauge group is U(l). We shall refer to the 
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internal degrees of freedom as colour degrees of freedom, and the associated 
gauge charges as colour charges. 

In what follows, quantities that possess colour are represented with Frak- 
tur characters; colour indices, when needed, will be given by lower case Fraktur 
characters, and the adjoint representation will be assumed. The gauge covariant 
derivative (both gauge covariant and covariant on the manifold A4) is given by 
(®At)% = V M <5% + f a b C 2l M c wnere 21// is the connection and f a t, c are the struc- 
ture constants of the group. The curvature of this gauge covariant derivative 
operator is the field tensor: 3^,/ = 2V[ M 2l„] a + f a b c 2l//2t,/. This field tensor is 
covariant under gauge transformations of the form 21// fy] = 2l//[0] + (£) M ) a by fa . 
We partially restrict the gauge freedom of the potential so that the components 
in an orthonorma! frame are finite everywhere in M. (except, perhaps, at truly 
pathological points such as curvature singularities). 

On a space-like hypersurface, S, it is possible to decompose the above quan- 
tities. The electric field as seen by an observer who is stationary with respect 
to the foliation is given by (£j a = h^^ u a u u . The magnetic field is given by 
23/ = — Tjjij i^v f " T -Sp < j a u v . The gauge covariant derivative that is compatible 
with the metric is (Ui) a b- 



3.1 The Dilaton- Yang-Mills Action and Variations 

The action is 

Sdvm= / ^v^J/vMwrj,;). (3.i) 

J M 4 

The colour indices are raised and lowered by the Killing metric g a b = ^fSafcb- 
The function, / YM (S f ), is a function of the dilaton alone, and contains no deriva- 
tives of the dilaton. 



Variation of the action of (3.1) yields source terms for the Einstein-Hubert 
field equations as well as for the dilaton field equation. In addition, variation 
with respect to the gauge potential 21// gives a source- free field equation for 
the gauge fields. Note that the variation, on the boundary, of the potential is 
the same as the variation of the potential projected onto the boundary. The 
induced variation in the Dilaton- Yang-Mills sector of the action is: 

SS DYM = f d 4 x^(-±(T DVM )^6g^ -±T DVM 5V + 

J M 

+ [ d 3 x(n DVM y a 5vu a - [ d 3 x(p DYM y a s^ a (3.2) 



where 

X% = (S) M ) 6 a (/ YM (W I/ 6 ). (3.3) 

This yields the Dilaton- Yang-Mills vacuum equation of motion, X v a = when 
the boundary terms vanish provided that there are no Yang- Mills source terms. 
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The Dilation- Yang-Mills stress energy and dilaton source are: 

*a) (3.4) 

and 

respectively. 

In addition we have momenta conjugate to the variation of the gauge fields 
on the boundaries T and S. These are: 

(n OYM )\ = V^U«m%^\n v (3.6) 

and 

(P^m)\ = -Vh = -VhM$)<? a (3.7) 

respectively. 

On the boundary T, we can decompose the variation of %i a into pieces 
normal and tangential to the foliation. Let 2J a = u 1 ^ and W a a = o l a *&i a . 
Then, 

(521/ = iV-VOKW) - 22J Q a <5Y a ) + a 4 a (52IT a a . (3.8) 

A similar decomposition of (7T DYM )a leads us to define a surface Yang-Mills 
charge density, 

(Qdym)o = v^/ym(^H<P„ (3-9) 
and a surface Yang-Mills momentum density 

(c7dym)<j = (Qdym ) 2Ha a . (3.10) 
Also, define a surface Yang-Mills current, 

V = v^/YM(^)e a TJ n^j a . (3.11) 
Then, the variation of 5 DY m on T is given by: 

SS DYM \ T = j &c ((J DYM ) a SN a - (Q DYM ) a 6(NV3 a ) + NT a 6W a a ). (3.12) 

We will consider the interpretation of the surface charge density in the following. 

3.2 Conserved Yang-Mills Charges 

In the presence of a source, 3^ a , the equations of motion for the Yang-Mills field 
(with dilaton coupling) are: 

(a M ) 6 a(/v M W i, 'b)=3V (3-13) 
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It can be seen that the source respects the identity (S AI ) b a 3 M b = 0. In an 
Abelian gauge theory (such as electro-magnetism) this quantity is gauge in- 
variant and can be used to define a conserved charge. However, in a general 
Yang-Mills theory, the identity is gauge-covariant, and the separation of the 
colour contained in the charge and the colour contained in the field it produces 
depends on the gauge choice. 

Yet, it is still possible to construct a conserved colour charge if we require the 
solution to the field equations to have certain properties (2^, |2j| . Suppose that 
the solution possesses a gauge Killing scalar, 6°, that is, a Lie-algebra- valued 
scalar field on M. that is covariantly constant: (£> M ) a b6 6 = 0. Then the quantity 
4°5 M o is gauge invariant and divergenceless: V AI (fi 3 Al a ) = 0. We can then define 
a charge, 

Qdvm[«] = - J^d 3 xVh Ufl t a T a (3.14) 

that is conserved provided that the source fCPo vanishes in the vicinity of T, 
that is, the charge of equation (3.14) has the same value regardless of the volume 
£ chosen for the integration. 

Recall now the equation of motion ( [3.13 ). Contracting both sides with the 
gauge Killing vector, one can take t a through the gauge covariant derivative to 
produce a gauge scalar as its argument. Further contracting both sides by 
and recalling that is proportional to a gradient, = —Nd^t, we can write 
n) b *{J yhWFs) as (^^(iW^W"^), due to the antisymmetry of 
the field tensor a in fi and v. Integrating over the space-like hypersurface 
S, we obtain 

QdymM=/ d 2 x(Q DYM ) a t a . (3.15) 

JdZ 

This alternate expression for the conserved gauge charge gives us our interpre- 
tation of (Qdym)cx as a surface gauge charge density. 



3.3 Canonical Form of the Dilaton- Yang-Mills Action 

Finally, we turn to the task of writing the Dilaton- Yang-Mills action of equation 
(3.1) in canonical form. A straight-forward calculation shows that £ t 2l ? a = 
-&j°t J - (5») a i-2li b where 2l t a = fSt M a . Using the decomposition of f into 
and N 1 *, we have 2t f a = iVQJ a - Y a 2rr a a . Then, we can show that 

(-^dym) a^t^i = ^(-Pdym) a^T (-^dym) cl$ij 

+% a (d,) b a (P DYM )\ - V,(2l 4 a (P DYM )\). (3.16) 

We also use the decomposition: 

jO,*' = - <£W)- (3.17) 
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The Hamiltonian density is given by H OYm = (P D YM) i a£t2lj a — L Dy - M where L D ym 
is the Lagrangian density (the integrand of equation (]3 . l[) ) . The Hamiltonian 
is, then, 



d X HnYM 

/ d 3 x(H DYM DYM 
J £ 

DYM 

) a X a N-(J DVM ) a N a ). (3.18) 



where 

W DYM - -(Pdym)\^ + \sfh / Y M(^)(» T a ^ a " (3.19) 

is the contribution to the Hamiltonian constraint from the Dilaton- Yang-Mills 
sector, 

dym)j — (-^dym) a'Sij (3.20) 

is the contribution to the momentum constraint from the Dilaton- Yang-Mills 
sector, and 

(<?DYm)o = o(^DYM)'b (3.21) 

is the Gauss constraint for the Yang-Mills field. 
The action in canonical form is simply: 



^DYM / ^ % ((-PdYm) a^t^l -^D 

J M 

= j dtUd 3 x((p DYM y a mi 



—HdymN — {TL m - M )iN l — (0DYM)a2l* a ) 

+ / d 2 x(~(Q DVM ) a WN+{J DVM ) a N a )). (3.22) 

JOE J 

4 Axion Sector 

In this section we consider the coupling of an axion field to the Yang-Mills field 
strength. Such couplings provide an interesting counterexample |15| to Schiff's 
conjecturco |24|] and, in the case of an Abelian gauge theory, imply interesting 
new tests of the equivalence principle |25, 26 . 



Many of the derivations in this section are similar to those in the previous 
section on the Yang-Mills field. Define the dual to the Yang-Mills field tensor 
by (*3V° = ^ pa d P a a . Note that there is the identity: (£> J" a (*3)^ b = 0. 



2 Schiff's conjecture states that any self-consistent theory of gravity that obeys the weak 
equivalence principle necessarily obeys the Einstein Equivalence principle 
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4.1 The Axion-Yang-Mills Action and its Variation 

We take the action for the Axionic sector to be 

S AYM = [ d 4 a;V^(^ YM (0)(^)^ Q ^ Q + ^ KE (0)(V0) 2 +^ PE (0)) (4.1) 

J M 

where 9 is the axion field and i? YM (0), i? KB (0), and i? PE (0) are functions of the 
axion (but not to its derivatives) that are the couplings to the Yang-Mills fields, 
the Kinetic energy, and the Potential energy respectively. 

Varying the action with respect to the geometry, the gauge field, and the 
axion yields 

SS AYM = f A V^d(-l(T AVM )^ v Sg^ + 2T a JV + 
Jm 

+ d 3 x ((n AYM y a s% a + n axi 5e) 

- d 3 x ({P AYM y a 6%j a + P axi S0) (4.2) 

where the stress-energy-momentum contribution from this sector is: 

(T AYM )„„ = -2t9 ke (0)(V m 0)(V„0) +ff^7? KE (0)(V0) 2 +g^ PE (6). (4.3) 

The equation of motion for the gauge field, when boundary variations vanish, 
is 2)% = (vacuum case) where: 

?Ta = W.V^vmW ( 4 - 4 ) 
while the equation of motion for the axion itself is — with: 

w _ l d^je) B muo) 2 2 dtf PB (g) 

" axi -4^T~ ( * y) °^ + _ d0~ « W) " 2V (45) 

Similarly, there are momenta conjugate to the gauge field configurations on T 
and E, 

(n AYM )\ = ^^m{0W%{^T\ (4.6) 

and 

(Paym)\ = -V^YM^K^^a (4.7) 

respectively. Also, there are momenta on T and S conjugate to the axionic field 
configurations: 

n axi = 2^$ KE (9)n^d^8 and P axi = -2v^tM0)u"0 M ( 4 - 8 ) 
respectively. 
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Recall that the variation of the gauge field on the boundary T can be de- 
composed as in equation ( |3.8[ ). We perform a similar decomposition of the 
momentum (-Z7 AYM )V Define a surface axion- Yang-Mills charge density, 

(Qaym)o = V^^vm(^)%®*o (4.9) 

a surface axion- Yang-Mills momentum density, 

) a 2U a ° (4.10) 

and a surface axion- Yang-Mills current, 

*a B = V^$YM(0W 3 n,£j a . (4.11) 
Also define the scalar density 

A = N- 1 n &xi . (4.12) 
Then, the T portion of the variation of S AYM is decomposed into: 

\t = J d 3 x ((J AYM ) a 6N a - (Q AYM ) a 6(NV3 a ) + N(M. a a 6W a a + A59)). 

(4-13) 

In analogy with the discussion of conserved charges in the previous section, we 
can interpret the quantity (Q A ym)o to be some sort of magnetic charge surface 
density that will yield a magnetic charge: 

Q DFF [t] = / d 2 x(Q AYM ) a t a (4.14) 

JdE 

when some gauge Killing scalar, 6 a , is present. 



4.2 Canonical Form of the Axion- Yang-Mills Action 



The Hamiltonian density of the Axion- Yang-Mills sector is given by H AYM = 
(PaymY a£t%k a + -Paxi£*# — L AYM . Denote by 9 the qu antity u^d^d. The first 
term in the Hamiltonian density looks just like equation ( 3.16 ) but wit h (Pa ym)^ 
replacing (P DYM Ya everywhere. The second term is just like equation ( 2.29 ), but 
here we replace Pdii and with P ax ; and 6. In addition, we can decompose the 
terms in the Lag icingiaii density, L^ymj 

which is just the integrand of equation 

(O). Note that: 



(4.15) 



Furthermore, the kinetic term of the axion can be decomposed just as the kinetic 
term of the dilaton was decomposed in equation ( [2.31 ). Thus, the Hamiltonian 
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is: 



d 3 xH k 



d 3 x(H AYM N+{H AYM 

+ I d 3 x ((Q AYM ) a X a N - (J AYM ) a N a ). (4.16) 
Js 



where 



n AYM = -(p AYM )\£, a + p a J-Vh(\$ YM {e)w a £ l a 

+$ KE (6)((V6) 2 - 9 2 ) + tf PE (0)) (4.17) 
is the Axion- Yang-Mills sector contribution to the Hamiltonian constraint, 

(Haym)j - -(P A ym)\ V + PaxlVjfl (4.18) 

is the Axion- Yang-Mills sector contribution to the momentum constraint, and 

(£a Y m)c = &)\(P A ™)\ (4.19) 
is the Axion- Yang- Mills sector contribution to the Yang-Mills gauss constraint. 



The action of equation (4.1) may be written in canonical form as 

Saym = / d X ((-Paym) a^t^i + Paxi^t^ — H A ym) 
J M 

= JdtU^ d 3 x ((P Ani )\m, a + P axi £ t 

—'HaymN — (^aym).^' — (G A YM)a%t a ) 

+ J d 3 x {-{Q AYM ) a X a N + (J AYM ) a N a )\ (4.20) 

5 Dilaton Four-Form Sector 

It is possible to treat a cosmological constant as a constant of motion arising 
from a four- form field rather than as a fundamental constant. Doing so allows 
one to consider space-times in which the cosmological constant takes different 
values. In this section, we will consider such a four-form field, but here we will 
also allow possible couplings to the dilaton. 

The four-form field strength will be given by A^ vpiJ — AW ^A vpa ^ where 
A\^ u is a three-form potential. The field strength is invariant under gauge 
transformations of the potential of the form A A/Jl/ [x] = ^4a/^[0] + 3!V[ A x M i/] 
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where x^v is an arbitrary two-form. We partially restrict this gauge invariance 
in requiring the components in an orthonormal basis be finite everywhere in M. . 
It will be useful to introduce an 'electric' field, Emj = h^h? hjA llup!J u' y . One 
can show that A^P a A puprJ = -AE h ^E m . 



5.1 The Dilaton- Four- Form Action 

Let /ppO?) be a function of the dilaton (that contains no derivatives of the 
dilaton) that couples to the four-form field as follows: 



$DFF — 



d'x^-g-^-f^A^Ax^p. (5.1) 

M 1 ' 4! 



We could replace the four-form field strength A\ pvp with its dual *A. The result 
would be a cosmological constant term with dilaton couplings. We will not do 
so here as it is more useful to consider variations of the three-form potential. 
We will discuss the relationship between the four-form field strength and the 
cosmological constant below. 

Under variations of the three-form potential, the geometry, and the dilaton 
field configurations, the induced variation in the action of equation (5.1) is: 

J M 

+ J d 3 x(n OFF )^ k SA ljk ~ J d 3 x(P DFF ) m U mj . (5.2) 

Note that the projection of the variation of the potential onto the boundary 
elements is the same as the variation of the projection of the potential onto the 
boundary elements. Here, the stress-energy-momentum of the four-form field is 
given by: 

(T DPP ) P „ = ^JW(^) (I g^A a ^ s A aPlS - A^A vaPl ) (5.3) 
and the dilaton source is: 

1 dfm(W) vpa . . 

The equations of motion for the four-form fields are: 

(~ DFF ) A ^ = -V a (f FF (<P)A aX ^). (5.5) 

In addition, the momentum conjugate to the three-form potential on the bound- 
ary T is given by: 

(n DFF r k = ^=7 U^)naA aX ^%& v (5.6) 
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while the momentum conjugate on the boundary S is: 

(P DPP )^~ = ~Vhf FF (^)u a A aX ^h h x ^hl = Vhf FF (<P)E h ^. (5.7) 

Since Aijk is proportional to the volume element e^/- on T, the quantity 
Vij = AijkU k will be a tensor on <9£, that is, Vij — hfh^Vab- With this in mind, 
we see that the variation of A^k can be decomposed as follows: 

5A ijk = SN-^Nu^ - 3u^a b k] 5V ab . (5.8) 

With the definition 

(Q DFF ) ab = 3^/ FF (^)n^ afc (5.9) 

for the four-form surface charge density, the component of the variation of the 
action S OFF on the boundary T can be written: 

SS DFF \ T = - J^d 3 x(Q DFF ) ab S(NV ab ). 
As usual, (Qdff)" 6 is related to a surface charge density. 



5.2 The Conserved Charge of the Four-Form Field 

Recall the equation of motion for the four-form field strength given by equation 
(|5~5|). In the presence of a source, J x,JlV , we have: 

- V Q (/ FF (f)yl QA ^) = J x ^ (5.10) 

and the identity Va J x w follows. Let the quantity be the volume element 
on the surface dS. Then, Vj(J A '" / e fIi/ ) = 0. If there is no source present in the 
vicinity of T, then we see that 

d 3 xy/hp (5.11) 

is a conserved charge of the four-form field: its value is independent of the choice 
of leaf SB Here, p = — 3e a b J^^u^ is the charge density. 

Using our definition of p, we can express it in terms of the four-form field 
strength by means of the equation of motion, ( 5.1 0| ) . Because the unit normal 



to S is proportional to a gradient (recall = —Nd^t), it can be pulled inside 
the derivative operator due to the antisymmetry of the field strength; in doing 



3 That is, given a two-form, e M „, the associated charge is independent of the foliation. 
However, the interpretation of this two-form as the volume element of the boundary of the 
leaves of the foliation is, of course, foliation dependent. 
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so, the derivative operator on M becomes the one on E. We find that yhp = 
3VhVi(f FF (&)E Tab e ab ) and thus, 

Q DFF [e] = / d 2 x(Q DFF ) ab e ab . (5.12) 
JdS 



Equation (5.12) justifies our interpretation of (<2 DFF ) as a surface density of 



four-form charge. 

5.3 Canonical Decomposition of the Dilaton- Four- Form 
Action 



To obtain the Hamiltonian density, -ff DFF = (P DPP ) J £t^ftjj— -L DPF , we compute: 

-3A m \7 h (P DFF ) m i + V h (3NV ab (P DFF ) hab l5.l3) 



(P DFF ) m J£ t A mj = N^-^yhf FF (>F)E h ^E m3 



Here, A^ t = A^t^. We have already seen that the Lagrangian density is purely 
electric: 

L DFF = -NVh J—f^{W)I^E m . (5.14) 
Thus, the Hamiltonian is: 

d 3 x H DF 



FF 

S 



NV ab (5A5) 



= f d 3 x (H DFF N + {g DFF fU m ) + f d 2 x (Q DFF ) a; 

J S JdS 

where the four-form contribution to the Hamiltonian constraint is: 

H UFF = -Vh^f FF (&)E m >E m (5.16) 
and the four-form Gauss constraint is: 

(Sbff) 53 = -3V ri (F DFF ) r ^. (5.17) 
The action in canonical form is: 

Stiff = / d^X (^(P DFF )^ £tAf^j — i? DFF ) 

Jm 

= j ' dt Qf d 3 x ((p DFF ) m3 £ t A h ,j - h dff n - (g DPP y>A w ) 

- [ d 2 x(Q DFF ) ab NV ab . (5.18) 

JdS 
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5.4 Relation to a Cosmological Constant 

Recall the Lagrangian density 

= -\V~9 UAm*Af . (5.19) 

In the latter form, it appears like a dilaton coupled to a cosmological constant 
A = — j(*A) 2 , except that, where a cosmological constant is fixed in the theory, 
*A is a scalar field. The stress-energy-momentum tensor can be evaluated and 
we find that (T DFF ) M „ = §/ PP (# ){*A) 2 

9yiv which is consistent with this interpre- 
tation. When the Gauss constraint equation holds, however, we find that the 
quantity f FF {^)Etnj must be equal to the tensor etuj up to some constant 2C. 
Thus, we find that / FF (!?)(*i) = 2C. The Lagrangian density is now written 
L OFF = — 1\/—g (l// FF (!f - ))C 2 . Thus, C 2 can be interpreted as a (negative) 
cosmological constant when f WF (^) is a constant. When / FF (<?) is not a con- 
stant, L DFF is essentially a contribution to the dilaton potential proportional to 
C 2 times the reciprocal of the coupling of the four-form field strength with the 
dilaton. Although C is a constant (rather than a scalar field), it is a constant 
of integration rather than an imposed theoretical constant. 

When the Gauss constraint holds, we can obtain an explicit expression for 
the charge Q DFF [e]. We see that (Q DFF ) ah = 6Cy/ae ab , so the charge Q DFF [e] is 
just 12C times the area of the surface dS whose volume element is e a b- 



6 Dilaton Three-Form Sector 

In low energy string theory, one encounters a term in the action that is a three- 
form field strength squared with a coupling to a dilaton. Here, we analyze 
such a gauge field with a somewhat more general coupling. The three form field 
strength will be written as = 3! V^A^i where A^ v is a two-form potential 

field. As always, there are gauge transformations of this potential, A p ^[x] = 
j4/h>[0] + 2! Vf^XHi under which the three-form field strength is invariant. As 
usual, we can use this gauge freedom to guarantee that the components of the 
potential in the orthonormal frame are finite everywhere in M.. It is useful to 
decompose the three-form field strength into 'electric' and 'magnetic' pieces on 
a space-like hypersurface, Ejj — h^hjH\ l _ Ll/ u' J and B = —{S\)~ 1 e x ^ a H\ l _ Ll/ u a 
respectively. It can then be shown that H X ^H\^ U = 3B 2 - QE^E^. 

6.1 The Dilaton-Three-Form Action 

As usual, the dilaton three-form action will be taken as the three-form field 
strength squared coupled to some function of the dilaton. Let this function be 
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given by f TF (\P), and make the usual assumptions about its form, i.e., that it 
contain only the dilaton, but no derivatives of the dilaton. Explicitly, we write 
the action: 

S DTF = [ d A x^-^—f^)H x ^H XlllJ . (6.1) 

J M 1 ' 

If we were to perform a duality rotation of the three-form field strength, we 
would be left with a kinetic energy for a scalar field. However, we are not 
interested in this case here. (In a sense, we have already considered a scalar 
field in the context of the dilaton in which there is no coupling to the Ricci 
scalar, though we have not considered a dilaton coupled to a scalar field.) 

The action of equation (|6.l|) can be subjected to variations in the geometry, 
the dilaton, and the two-form potential. The response in the action is given by: 

J M 

+ d \r (/7 DTP )tf Aij - J d 3 x (P DTP )^-. (6.2) 

As usual, the projection onto the boundary elements of the variation of the 
potential are the same as the variation of the projection of the potential onto 
the boundary elements. When the variation of the two-form potential on the 
boundary is held fixed, there is a contribution to the stress-energy-momentum 
tensor: 

(T DTF )^ = / TF (f)((2 • 3!)- V# Q/37 #«/37 - \H^H uaP ) (6.3) 
as well as a dilaton source: 

T mF = - l -^l H ^H x , v . (6.4) 

In addition, the equation of motion for the three-form field, in the absence of a 
source, is the vanishing of 

{S^y v = -Vx{f™{V)H x » v ). (6.5) 

The momenta conjugate to the two-form potential configuration on the bound- 
aries T and S are, respectively, 

(iWr = ^h.mnxH x ^f^l (6.6) 

and _ 

(■PdtfF = -VhUWuxH^fyhl = -Vhf TF (V)E^. (6.7) 

The potential Ay on the boundary T can be written in terms of the quan- 
tities V a = a^AijU^ and W a b = ^aPb-^i on ^ ne two-surface dS as 

Aij = 2u {i al V a + d%o% W ab . (6.8) 
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Note that since W ab is a two-form denned on the two-surface dE, we could 
express it as W a b = we ab where w is a scalar function on the two surface. 
Similarly, the variations of A y - can be decomposed: 

<L% = 2u [l a] ] N- 1 (S(V a N) - W ab SN b ) + &ffi } 6W ab . (6.9) 

This leads to a like decomposition of the momentum conjugate to the T bound- 
ary two-form potential. We construct a surface three-form charge density, 

(QDT F ) a = 2V^/ TF (f r K£ Aia (6.10) 

which is a vector on OS. Note that (Q D tf) q is an extensive quantity. It can 
be used to define conserved charges as will be seen below. Another extensive 
quantity is the surface three-form momentum density, 

(iTdtf) b = (QnT F TW ab . (6.11) 

Also, there is a surface three-form current density, defined as 

(X DTF ) Qb = ^/ TF (^KiJ^ h . (6.12) 

Using these, we can write the T boundary contribution to the variation of the 
dilaton three-form action as: 

SS UTF \ T = J^d 3 x ((Jo TF ) a 5iV Q - {Q OTF ) a S(NV a ) + N(l DTF ) ab SW ab ). (6.13) 

Here, the first two terms in the integrand involve variations of intensive variables 
with extensive coefficients while the last term is a variation of an extensive 
variable with an intensive coefficient. 

6.2 Conserved Charges of the Three-Form Field 

Suppose that suitable matter is present such that there is a source term, J^ u , 
included in the equation of motion of the three-form field. Then we have from 
equation ( |6.5| ), 

-V X (U(W)H X ^) = J"". (6.14) 

However, we see from this equation that the source must be divergenceless, 
V A1 J /il/ = 0. Thus, we can construct conserved charges for it in the following 
way. Let — V M be an exact one-form. (0 is a scalar function.) Due to 
the antisymmetry of in \i and v, we have = /^V^J^" = V ^{J^ f v )- 
If we integrate this expression over M, we have contributions from the initial 
and final hypersurface as well as the boundary T. Supposing that the source 
vanishes in the vicinity of T, we find that the quantity 

Qdtf[/] = - / (Fxy/hu^U (6.15) 
J s 
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is invariant upon the surface, S, of evaluation. Thus it is a conserved charge. 



Using equation (3.14), we can re-express the conserved charge of equation 



( |6.15| ) in terms of the two-surface three- form charge densit y (Q DTF ) a . Recall 



that = — N'Vfit. Then, the integrand of equation (6.15) can be written in 
the form V p (f TF (^)E^ a f a ) . (Notice that we have restricted the gradient f a to 
be a one-form on the dual tangent space of dS. Natural choices for the scalar 
<j) are the coordinates of this hypersurface.) Therefore, we have: 

Qdtp[/] = / d 2 x(Q DTF ) a f a (6.16) 

JdE 

where f a = d a <t> and <f> is an arbitrary scalar function on OS. 

6.3 Canonical Form of the Dilaton-Three-Form Action 

The canonical decomposition of the action >!5 DTF is obtained through the con- 
struction of the Hamiltonian density H DTF = (P DTF y^ £ t Ajj — L DTF . A straight- 
forward calculation of the first term yields: 

(p DTF yi£ t A- t3 = ±N(p DTF y3E- l3 + ±N h H m3 (p DTF yi 

-2A tJ V,(P DTF y3 + 2V,(A tJ (P DTF )^) (6.17) 

where A t j — t^A^j will act as a Lagrange multiplier. Also, the second term in 
the Hamiltonian density can be written as: 

l dtf = N^h / TF (<f)(is 2 - ^ /•:"/•;, j. (6.18) 

With these expressions, we can write the Hamiltonian: 

TF 

= / d 3 x (h dtf n + (n DTF yN* + A tJ {g DTF y) 

Js 

+ f d 2 x((Q DTF ) a NV a -(J DTF ) a N a ). (6.19) 

JdE 

Here the contribution from the dilaton three-form sector to the Hamiltonian 
constraint is given by: 

H DTF = -\VhB 2 (6.20) 
the contribution to the momentum constraint is: 

(Wdtf)r = \H Ulj [P^Y 3 (6.21) 
and the Gauss constraint for the three- form field is: 

(g mF y = 2V,(P DTF f? (6.22) 



E 
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Also, the action of equation (6.1) can be re-expressed in the canonical form. It 
is: 



*Sdtf — I d X ((-£dtf) ^ ^t-A-ij ^dtf) 

Jm 

= J dt(^J d 3 x ({P DTF yi£ t A- tJ 



— HdtfN - (Hdtf)iN 1 — (5 D TF)'Ai) 

+ / d 2 x ((J DTF ) a N a - {Q DTF ) a NV a )) . (6.23) 

JdS J 

7 Statistical mechanics 

Here we address, at a formal level, the construction of statistics of the quantum 
mechanical theory based on the actions we have constructed. The mechanical 
theory is entirely described in terms of path integrals involving our actions. 
These yield quantum mechanical density matrices. The statistics of primary 
interest are the partition functions, which can be thought of as functional in- 
tegrals over the density matrices with some periodic identification (with the 
period related to the temperature). 

7.1 The General Action and Canonical Ensembles 

We will consider a general form of action that possesses all of the features of the 
actions we have considered thus far. Let A^ A denote the gauge field potential, 
where the upper-case Latin index labels the various gauge fields that may be 
present. The field strength associated with the potential is given by A , and 
the gauge field action functional, S GF , involves the square of the field strength 
along with a possible coupling to the dilaton, / GF (!f r ). Note that for the three 
and four form field strength gauge fields, the field index will actually be a set 
of space-time indices, while for the Yang-Mills fields, it is a colour index. The 
total action of the theory is given by Stotai = Sehd + S GF + S . Variations of 
this action yield the usual equation of motion terms as well as boundary terms. 
We will primarily be concerned with the latter. The T portion of the variation 
of the total action will take the form: 

<5S total | r = J T d 3 x (-ESN + J a SN a - Q A S(NV A ) 

+N(S ab 5a ab + y5V +l a A SW a A )). (7.1) 

Here, Q A are the surface charge densities of the gauge fields, T a A are the surface 
currents of the gauge fields, and J a = (Jehd)o + QaW u a is the total surface 
momentum density. Also, V A are the projections of the T boundary potentials 
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of the gauge fields along u l , while W a A are the portions of these fields restricted 
to dS. Note that all of the gauge field actions we have considered can be 
expressed in this way. The action can be written in canonical form as: 



Stotal = J dt (J^ d A X (pVlthv + P da £ t # + P\£tA A 

-UN -H-iN* -G A A t A ) 
+ [ d 2 x(-£N + J a N a - Q A NV A )\ (7.2) 

JdS J 

Here, A t A are the projections of the gauge potentials along the time vector. H 
is the total Hamiltonian constraint, Ti. a is the total momentum constraint, and 
Q A are the Gauss constraints of the gauge fields. 



The first three terms in the integrand of equation (7.1) involve variation 
of an intensive variable with an extensive coefficient. However, the remaining 
terms involve variations of extensive variables with intensive coefficients. Define 
the micro-canonical action to be 

Smicro = Stotai + / d 2 x (£ N - J a N a + Q A NV A ) . (7.3) 

JdS 

Since this action differs from the original by boundary terms alone, the equa- 
tions of motion are unaffected. However, such a transformation changes the T 
boundary component of the variation of the action. We now have: 

SS micro \ T = J^d 3 xN(6£-Lo a 6J a + V A 6Q A 

+ \S ab 5<y ab + y6& + l a A SW a A ) (7.4) 

where we define Nuj a = N a . We will interpret io a as the angular velocity of 
observers co-moving with the foliation on T. Each term in the integrand of 



equation (7.4) is of the form of a variation of an extensive variable with an 
intensive coefficient. We see that the boundary conditions imposed in order to 
obtain the equations of motion are those in which the extensive variables are held 
constant. It is for this reason that the term 'micro-canonical' has been adopted. 
Also, note that the Hamiltonian for the micro-canonical action contains only 
the constraints (i.e., there are no additional boundary terms); so the canonical 
form of the micro-canonical action has no dE integral. The importance of this 
observation is that the micro-canonical action vanishes when stationarity of all 
the fields and the constraint equations are imposed. 

Grand-canonical boundary conditions involve the fixation of all intensive 
variables on the boundary. In analogy with the procedure above, we define a 
grand-canonical action by: 

5 grand = S to tai - J^d?xN(\S ab a ab + y# +J a A W a A ). (7.5) 
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Under variations of this action, the T boundary component becomes: 

JSgrandk = - J &c (ESN - J a 5{Nuj a ) + Q A S{NV A ) 

+a ab 5(NS ah /2) + V5(Ny) + W a A S(Nl a A )) . (7.6) 

so it is indeed the intensive boundary variables that must be fixed in order to 
obtain the equations of motion. In canonical form, the T contribution to the 
grand-canonical action is given by: 

Sgrandlr = - j d 3 x N (y - J a uo a + Q A V A + \ tr(S) + W + I a A W a A ). (7.7) 

It will be useful to have a covariant form of the micro-canonical action. 
To obtain this, we recall the definition of the micro-canonical action, equation 
(f7Tg|), as well as the initial covariant form of the total action, equation (2.1) plus 
matter terms. Since the difference between the micro-canonical action and the 
total action amounts to just a difference on the boundary T, we will just attend 
to this. Thus: 

Smicrolr = J <fic (-2V^7/E H (f ) tr(0) - N£ + N a J a - NV A Q A ). (7.8) 

To cast this in a covariant form, we employ the definitions of £ , J a and Q A , as 
well as the decomposition given in equation (2.10). Then equation (7.8) can be 
re-written in the covariant form: 



Smicrok = / d 3 x^(f EH (<P)t fl 6>^d„t-2n»dJ EH (<P) 



,3 ~ 

7 

-f G¥ m{n^ A d v tWA A )). (7.9) 



7.2 'Euclidean' Notation 



Until now, the metrics we have considered have been of Lorentzian signature. 
Such metrics can be re-expressed in a 'Euclidean' notation by redefining the 
(real) metric functions in terms of new (complex) ones. Our prescription is the 
following: We first rewrite the volume elements as ^fg — iy/~g and ^Fy = iy/— 7 
and redefine the lapse and shift functions: N = iN and N l = iN l . The former 
leads to a new unit normal = — iu M which satisfies u ■ u = +1. In addition to 
the lapse and the shift, we adopt a new notation for all the Lagrange multipliers 
in the Hamiltonian. Thus, A t A = \A t A . Thus, the Hamiltonian as a functional of 
the Lagrange multipliers is B.[N, N\ A t A ] = iH[JV, N l , A t A }. Finally, we define 
a new action functional I[g, <?■, A] = —iS[g, \P, A] such that the phase in the path 
integral is written exp(— I). 
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An important feature of this prescription is the following. Suppose that a 
complex metric g c is obtained from a real, Lorentzian metric g L via the trans- 
formation t — > —it. Then, the Lagrange multipliers, N, iV l , and A t A become 
imaginary, or, equivalently, N, N l , and A t A become real. However, the canon- 
ical data, that is p 1 ^, h^j, P^n, \P, P % A and A^ A all are invariant. Therefore, 
the extensive variables are all invariant under the Wick rotation since these 
variables are constructed out of the canonical data. In particular, the values of 
the extensive variables of the complex metric that extremize the path integral 
are the same as the values of these variables on the corresponding Lorentzian 
metric. 

We can write the micro-canonical and the grand-canonical actions using the 
Euclidean notation in the canonical form. They are: 

/micro = J dt J d 3 X (-ip«£ t % - iP\ £ t A, A - LP dil £ t <f 

+HN + H T N l + G A A t A ) (7.10) 



and 



/grand = /micro + dt d?X N (£ - J a UJ a + Q A V A + 

J JdS 

\tr{S)+y* + l a A W a A ). (7.11) 

In addition, we can write out the variations of these actions. For example, the 
T boundary contribution of the variation of the Euclidean form of the micro- 



canonical action can be constructed from equation (7.4). It is: 

<5/ ml cro|r = - [ dt [ d 2 xN(S£ ~uj a SJ a + V A SQ A 

J JdS 

+\S ab 8<j ab + y5V+I a JW A ). (7.12) 



7.3 Functional Integrals 

It is well known that the classical equations of motion are obtained from an 
action /, by setting variations of this action to zero. As argued earlier, the vari- 
ations are of a restricted class which fix certain quantities on the boundary, and 
the quantities that must be fixed are determined by the choice of ensemble. For 
example, it is the extensive boundary variables that must be fixed in obtaining 
the equations of motion from the micro-canonical action. 

This method of obtaining classical equations of motion is best understood 
as a (classical) limit of the quantum mechanical density matrix obtained from 
a path integral. Formally, we can write this density matrix as: 

q= J dfolMle-'b'*^ (7-13) 
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where the integral is taken over all possible field configurations of the metric, 
the dilaton, and the matter fields, from an initial and a final space-like surface. 
The field configurations are fixed on the initial and the final space-like surfaces, 
^initial and £s na i; so the density matrix is a functional of these quantities. 
In addition there will be quantities on T joining the boundaries of these two 
space-like surfaces that must also be fixed. These quantities depend on the 
ensemble chosen for the action. The density matrix will be a functional of 
these variables as well. The density matrix gives us all the quantum mechanical 
information we need to solve most problems. Furthermore, the classical limit 
can be understood in the following way. Given an initial configuration, there is a 
final configuration that is most 'likely' in a quantum-mechanical sense: it is the 
one for which the density matrix has the largest value. In the classical limit, the 
action is always much larger than unity (in units of Planck's constant). Thus, 
the final configuration that extremizes the density matrix is the one for which 
51 = 0, that is, the one for which the classical field equations hold. (Note that 
/ is a complex quantity — we have not yet made a Wick rotation — so the last 
statement follows from stationary phase arguments.) 

Much information about the system can be obtained by taking statistics of 
the density matrix. The most important of these is known as the partition 
function. It is obtained by 'tracing over' the initial and final states. In terms 
of path integrals, this can be realized by identifying the initial and the final 
configurations, and integrating over all possible configurations. The partition 
function then depends upon the information fixed on the boundary T alone. 
Note that, in identifying the initial and the final space-like surfaces, we have 
effectively changed the topology of M. from £ x X to £ x S 1 . 

As an example, consider the micro-canonical ensemble. The objects that are 
held fixed when this action is varied are the metrics and field configurations on 
the initial and final surfaces as well as the extensive variables on the boundary 
T. Thus, the functional dependence of the density matrix is: 

d[g,V,A]e- I ™^> 9 ' A \ (7.14) 



The subscripts 'i' and 'f refer to the initial and final space-like hypersurfaces, 
-S'initiai and £finai, respectively and ip is the dilaton configuration on T. The 
partition function is also known as the density of states. It is given by: 

u[8,J,(r,il>,Q,W\ = J d[h,<P s ,A s ] Q[h,# s> A s ;h,& s ,A s ;£,J,<j,ip, Q 7 W}. 

(7.15) 

Note that we have explicitly indicated the periodic identification of the fields 
on the initial and final space-like boundaries in the density matrix. 

In a similar manner, we can construct a grand-canonical density matrix by 



making a Laplace transform of equation (7.14). This density matrix would be a 



Quasilocal Thermodynamics. 



33 



functional of (3, ui a , S a b, y, V A , and I A on the boundary T. Here, the reciprocal 
temperature, /3, is defined by: 

p = jdtN\ as (7.16) 

where we have used § as a reminder that the initial and final space-like hypcr- 
surfaces have been periodically identified; (3 is the gauge invariant measure of 
this period. The grand-canonical partition function is then u>, S, y, V,T]. 



8 Black Hole Thermodynamics 

The thermodynamics of systems containing black holes can be obtained using 
the statistical techniques of the previous section. In particular, we wish to obtain 
a thermodynamic 'first law' which relates the entropy change of a system with 
changes in extensive quantities of the system. Clearly, such a result should be 
obtained from the micro-canonical action. The entropy is normally defined as 
the logarithm of the density of states. (We recall that the latter is just the 
partition function for the micro-canonical ensemble.) We will only consider 
black hole solutions that are stationary in the sense that the Lie derivative of 
all fields (including the metric, dilaton, and matter fields) vanish. Also, we 
continue to use the 'generalized gauge fields' introduced in the previous section. 
Any of the gauge fields discussed in this paper will have a similar form. 

We will only construct the 'zeroth order' contribution to the entropy from 
the path integral for the density of states. In this approximation, the density of 
states is written: 

v[£, J, cr, ip, Q, W] w c~ /micro (8.1) 

where the micro-canonical action is evaluated at its (complex) extremum value. 
Clearly, the entropy is just: 

^micro | cl • 

(8.2) 

Thus, to evaluate the entropy is to evaluate the micro-canonical action at its 
(complex) extremum value. However, the presence of an event horizon intro- 
duces some interpretational problems. These we address in the rest of this 
section. 



8.1 Regularizations at the Event Horizon 

As usual, we turn to the Euclidean form of the action. The system will be defined 
as the interior to some closed two surface, B outC r, on which the micro-canonical 
boundary data are specified. This interior will include an event horizon; the 
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foliation will become degenerate on this horizon. In addition to the usual fo- 
liation in the time-like direction, we will suppose that it is possible to foliate 
from the event horizon to the outer boundary. After periodic identification, the 
complex manifold has the form of a conexS 12 , where the azimuthal direction 
of the cone is the time-like direction and the 'radial' direction corresponds to 
the space-like vector along which the space-like foliation was constructed. The 
outer edge of the cone is the two-surface B outer , and the 'point' of the cone is 
where the time-like foliation becomes degenerate: the event horizon. 

It is common to impose regularity conditions on the cone to cause it to 
become a disk. Near the event horizon, we assume that it is possible to write 
the metric in the Euclidean form: 

dsl ~ N 2 dt 2 + M 2 dr 2 + cj ab dx a dx b . (8.3) 

Here, r is the coordinate of foliation in the space-like direction, and M is the 
analog of the lapse function defined as n M = Md^r. The conical singularity is 



not present in the metric of equation (8.3) if the 'circumference' of the circles of 



constant t have the value of 2tt times the proper radius near the event horizon. 
Let At be the period of the identification of the coordinate t. Then the circum- 
ference of the circles of constant t is given by NAt, and the regularity condition 
is At(n l b\N) = 2ir at the event horizon. Choose the period At to satisfy this 
condition. Then it takes the value At = 2tt/k h where = h^{b\N)(djN) 
(evaluated on the event horizon) is the surface gravity of the event horizon. 

8.2 Evaluation of the Micro- Canonical Action 

Although the event horizon is in no way distinctive to an infalling observer, 
nevertheless, for a system observer (i.e., one who lives on the system boundary 
S ou t e r) it represents a one-way membrane onto which information can approach, 
but from which nothing can come (at least classically). For such an observer, 
it would be inappropriate to attempt to calculate an action that included the 
event horizon. So let us remove an open set surrounding the event horizon 
infinitcsimally close. The remaining portion of the manifold now contains a new 
'boundary' However, this is not a system boundary on which thermodynamic 
data must be specified — it is just a tool that will represent the information lost 
in 'throwing away' the event horizon from the system. With this understanding, 
we shall denote this surface as B lnnQ:r . 

Let us recall the covariant form of the micro-canonical action. It is the 



usual action but with the 7^ u t r boundary term given by equation (7i)). There 
is no supplemented boundary term for the inner boundary. Using the usual 
techniques to perform a canonical decomposition of this action, and keeping 
track of the terms appearing on the new boundary B lnncri we find that 

J dt d 2 x(2^ f EH (V)n%N - 2sfd Nn%f Ea {&) 
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-N a J a + NQ A V A ). (8.4) 

in the 'Euclidean' notation. Here, /* is the thermodynamic action which differs 
from the micro-canonical action, /micro, in that the event horizon has been 



removed from the micro-canonical system. In obtaining equation (8.4), we have 
used the relation % = N^ 1 d i N. The action / m i cro was given in equation (7.10). 
Thus, we notice that the only boundary contribution to /* is on the boundary 

dinner- 

We wish to evaluate the action /* on the solution in which the field equa- 
tions hold — the extremal solution. The contribution from the action /micro of 



equation (7.10) vanishes: the terms involving Lie derivatives of the fields vanish 
due to the stationarity of the solution, and the remaining terms vanish because 
the constraint equations must hold on a classical solution. Thus the sole contri- 
bution to the action /» comes from the boundary contribution £>i nn cr of equation 
( |8.4| ) which arises due to the discarding of the event horizon from the system. 
In evaluating this term, we invoke the regularity conditions of the previous sub- 
section: N = 0, N a = 0, and J n l b\Ndt = —2ir. (The shift approaches zero as 
we approach the event horizon since we have chosen a zero-vorticity observer. 
The negative sign in the third equation arises since the normal vector is now 
inward directed, rather than outward as it was in the previous sub-section.) 
Thus we find that: 

J*| CL = -47r/ d 2 x^f EH (V). (8.5) 

JdS 

The entropy (neglecting higher-order quantum corrections) is just the negative 
of this quantity according to equation ( |B.2| ) (using /* rather than /micro)- The 
usual value for the entropy with no dilaton field — one quarter of the event hori- 
zon area — is recovered when we set / EH (^) = (167r) _1 . (This is the appropriate 
value for this coefficient for the usual Einstein-Hilbert action in units where 
Newton's constant is unity.) 



8.3 The First Law of Thermodynamics 

Since the inner boundary, Sinner, is not really a boundary, we will view /* as 
a functional of the same extensive variables on yB ou ter as /micro- The first law 



of thermodynamics is obtained by varying the entropy given by equation (B.2) 



with the aid of equation (7.12) which applies equally to /*. These variations are 
understood to be amongst those that preserve the classical equations of motion, 
so the boundary contribution is the only one present. (Recall that we have 
closed the manifold with respect to the initial and final hypersurfaces, so there 
are no longer any initial and final space-like boundaries.) Thus: 

5§= f d 2 xp(5£ ~u a 5J a + V A SQ A + ±S ab S<7 ab + ySV +l a A 5W a A ) (8.6) 



Quasilocal Thermodynamics. 



30 



is our formulation of the first law of thermodynamics. Recall that (3 is the 
reciprocal temperature of the system (equation (7.16)). Note that it is not 
necessarily constant over the system boundary, so we have left it within the 
integral. When the boundary is chosen to be an isothermal surface, then we can 
integrate the first term in the integrand of equation (S.6) to obta in the usual 
'/35E' term. (Recall that E is the quasilocal energy of equation ( 2. 16| ) . ) However, 
the angular velocity uj a will can not generally be a constant on the system 
boundary simultaneously, so the usual expression for the first law does not 
generally hold in the case of a finite system. In the case of non-abelian fields, if 
a surface can be chosen such that (3V A is constant and proportional to the gauge 
Killing scalar (if such exists), then the contribution of this term in th e first law 
is of the form pVSQ-p m where the latter quantity is defined in ( 3.14 ). We note 
that the formulation ( |8.6| ) of the first law differs from that considered previously 
for non-abelian gauge fields pjj ], in which certain asymptotic properties of the 
gauge and gravitational fields were assumed in order to define a colour charge — 
the resultant formulation of the first law is therefore valid only in this asymptotic 



region. Our formulation (S.6) of the first law recovers this result when the gauge 



and gravitational fields have the aforementioned falloff properties. 

It is useful to divide the variations of the metric on Z? ou t e r into a 'shape' 
preserving piece and a 'volume' (which is, of course, an area on a two-surface) 
preserving piece as follows: Sa a i, = ^ab^V^ + v^^afc, where s a b = °abl\f^ ■ If 
this is done, then the ^pS ab Sa a b term in the integrand of equation ( |3~6| ) can be 
re-written as (3(P5V + X ab S<;ab)- Here, V = ^fa can be thought of as a measure 
of the volume of the system (by which we really mean the area of B on t el ), and 
thus V = \ tr(S)/y/a is interpreted as the pressure on the system. The quantity 
A afc = ^y/aS ab is thermodynamically conjugate to shape changes of the system 
boundary. 



9 Concluding Remarks 

The formulation of the first law of thermodynamics as given in ( |8.6| ) generalizes 
previous formulations |L8| [2£| to include the most general couplings of gauge 
fields to dilatonic gravity in four spacetime dimensions that have at most two 
derivatives in any term in the action. We close by commenting on possible 
extensions of our work. 

General arguments from string theory suggest that the action considered in 
this paper receives corrections from terms that have more than two derivatives 
in the metric and matter fields, i.e., terms which are at least quadratic in the 
curvature and/or field strengths. Although the general form for the Noether 
charge for such terms has been evaluated p0| , the detailed manner in which 
such terms contribute to the first law of thermodynamics in the context of the 
quasilocal formalism considered in this paper remains to be worked out. 
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A related problem of interest concerns the role of topological fields. These 
are fields which couple to the connection but not to the metric. It has recently 
been shown that interesting black hole solutions exist for a model topological 
field theory in (2 + 1) dimensions p9| |. Their thermodynamical properties are 
considerably different from the usual case (l^, |3(| and are not fully understood. 
The generalization of the quasilocal formalism to such actions represents an 
interesting problem since the first derivatives of the metric will play a markedly 
different role in the boundary terms. 
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A Notation 

This is a summary of the notation used in this paper. The manifold M. is 
topologically S x X where S is a spacelike hypersurface and I is a timelike 
interval. Foliation of M along 1 allows us to define the lapse, N, and the shift, 
N^. The various manifolds we will consider, and some of the tensors defined on 
them are summarized in Table 1. We can construct the following densities on 

Table 1: Manifold Variables 



Manifold 





M 


r 


S 


dz 


Indicies 


{H,v,...} 




{M,J} 


{a,b} 


Normal Vector 






Up 


Ui 


Metric 




Hi 






Compatible Derivative 


v M 


Ai 


v, 




Intrinsic Curvature Scalar 


R[g\ 




R[h] 




Extrinsic Curvature 




®H 




k a b 


Geometric Momentum 











OS out of projections of 7r y : the surface energy density, £, the surface geometric 
momentum density, (J7ehd)<u the surface stress density S ab . 

We consider various forms of matter. These are summarized in Tables 2 
and 3. Table 2 contains the field variables used in the various sectors. The 

Table 2: Field Variables 



Field Strengths Potential Fields 

full electric magnetic full along u on dS 



Dilaton 




-V- 










Axion 




-9- 






-9- 




Yang-Mills 


1$ pv 


£i a 




V 






Four-Form 








A-\pv 


Vab 




Three-Form 


H\pv 




B 




V a 


w ab 


General Field 


F A 

± pu 






A A 


V A 


w a A 



momenta corresponding to these fields (and the related densities constructed on 
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dS) are given in Table 3. 

Table 3: Conjugate Momenta 



Conjugate Momenta Surface Density 





on E 


on T 


charge 


momentum 


current 


Dilaton 


Pdii 


n dil 




-y- 




Axion 


1 axi 






-A- 




Dilaton- Yang-Mills 


(-Pdym)^ 


(-^dym) Z ci 


(Qdym)o 


(uTdym ) a 


-J 


Axion- Yang-Mills 


(Paym) a 


(-^aym) a 


( Q AYM ) a 


(iTAYlja 


Jt a 


Dilaton- Four- Form 




(tf DPP p 


(Q DFF ) ab 






Dilaton-Three-Form 


(Pdtf) ?J 


(tf DTP )« 


(QDT F ) a 


(j7oT F )a 




General Field 


p\ 








ja 
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